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We derive the virial theorem appropriate to two-dimensional point vortices at statistical equi- 
librium in the microcanonical and canonical ensembles. In an unbounded domain, it relates the 
angular velocity to the angular momentum and the temperature. Our expression is valid for an 
arbitrary number of point vortices of possibly different species. In the single-species case, and in 
the mean field approximation, it reduces to the relation empirically obtained by J.H. Williamson [J. 
Plasma Physics 17, 85 (1977)]. 

PACS numbers: 47.10.A-,47.15.ki,47.32.C- 

I. INTRODUCTION 

In a seminal paper, Onsager [l| laid down the foundations of the statistical mechanics of vortices in two-dimensional 
hydrodynamics. He considered the point vortex gas as an idealization of more realistic vorticity fields [Tol ] and 
discovered that negative temperature states are possible for this system. At negative temperatures, corresponding to 
high energies, like-sign vortices have the tendency to cluster into "supervortices" similar to the large-scale vortices 
(e.g. Jupiter's great red spot) observed in the atmosphere of giant planets. The qualitative arguments of Onsager 
were developed more quantitatively in a mean field approximation by Joyce & Montgomery 0, 0|, Kida [|| and 
Pointin & Lundgren [H, 0] , and by Onsager himself in unpublished notes Q ■ The statistical theory predicts that the 
point vortex gas should relax towards an equilibrium state described by the Boltzmann distribution. Specifically, the 
equilibrium stream function is solution of a Boltzmann-Poisson equation. At positive temperatures, the Boltzmann- 
Poisson equation is similar to the one appearing in the theory of electrolytes in plasma physics (like-sign vortices 
"repel" each other) 0]. At negative temperatures, it is similar to the one appearing in the statistical mechanics of 
stellar systems (like-sign vortices "attract" each other) 0, [l(| • Many mathematical works (ITl - fl6| have shown how a 
proper thermodynamic limit could be rigorously defined for the point vortex gas (in the Onsager picture). It is shown 
that the mean field approximation becomes exact in the limit N — > +oo with 7 ~ 1/N, where N is the number of 
point vortices and 7 their individual circulation. The kinetic theory of the point vortex gas has been developed in 
Refs. [T3-H!- 

The statistical equilibrium state of a single species system of point vortices in an unbounded domain was studied 
by Williamson (2(| . He used a mean field approximation and numerically solved the Boltzmann-Poisson equation for 
an axisymmetric distribution of point vortices. He empirically obtained a relation between the angular velocity, the 
angular momentum, the number of vortices and the temperature but did not manage to derive it. The same relation 
was obtained independently by Lundgren & Pointin @ who derived it from the microcanonical distribution of the 
point vortex gas by using essentially the same trick as the one used by Salzberg & Prager [27], [H| to obtain the exact 
equation of state of a 2D plasma. On the other hand, Kiessling [231 derived the Williamson algebraic relation by 
a series of partial integrations for radial solutions of the Boltzmann-Poisson equation (he proved that the solutions 
are necessarily radially symmetric with decreasing density profile). He mentioned the similarity with a virial identity 
obtained by Corngold [30( but also stressed that this identity is valid only at positive temperatures and corresponds 
to different types of averages. Therefore, its application to point vortices at negative temperatures is not clear. 

In this paper, we provide a simple derivation of the virial theorem of point vortices at statistical equilibrium and 
show that, for an unbounded system, it is equivalent to the Williamson relation. Our derivation is more general than 
previous ones since it is valid for a multispecies gas of point vortices and goes beyond the mean field approximation 
(furthermore, it does not assume, or use, the fact that the flow is axisymmetric). Therefore, the relation that we 
obtain is "exact" . We also obtain a generalization of the virial theorem for an axisymmetric system of point vortices 
enclosed within a disk in the mean field approximation. However, in that case, the equation is not closed (as it involves 
the pressure on the boundary). Finally, we mention numerous analogies between 2D point vortices and 2D rotating 
self-gravitating systems. In particular, the statistical equilibrium state, the local equation of state, the condition of 
hydrostatic equilibrium (in the rotating frame) and the virial theorem have a similar form. The virial theorem plays 
a very important role in the theory of simple liquids [3l| and in astrophysics [32j . To our knowledge, the "virial of 
point vortices" (especially at negative temperatures) has not been introduced before. In view of its importance in 
other domains of physics, we think that it is important to devote a specific paper to this topic. 

The paper is organized as follows. In Sec. HH we recall the basic equations of the multi-species point vortex gas 
in an unbounded domain. In Sec. IIIIl we consider the statistical equilibrium state of point vortices in canonical 
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and microcanonical ensembles, and derive the Yvon-Born-Grecn (YBG) hierarchy and the "exact" virial theorem. 
In Sec. IIV[ we consider the mean field approximation, and derive the multi-species Boltzmann-Poisson equation 
(from the YBG hierarchy or from a maximum entropy principle) and the mean field virial theorem. We also consider 
some applications of the virial theorem for point vortices in an unbounded domain or in a disc (for axisymmetric 
distributions). The Appendices contain several additional results. In Appendix [XJ we define and compute the virial 
of point vortices. In Appendix [B] we provide alternative derivations of the virial theorem in an unbounded domain 
starting directly from the partition function or the density of states, or from the explicit analytical solution of the 
Boltzmann-Poisson equation (for a single species system). In Appendix [C] we derive the virial theorem in a bounded 
domain for purely logarithmic potentials (i.e. neglecting vortex images in the case of point vortices), and comment 
on the meaning of the "pressure of point vortices" . In Appendix [Dj we derive the two-body correlation function and 
the virial theorem for a neutral spatially homogeneous system of point vortices. Finally, in Appendix [El we discuss 
the "exact" values of the critical inverse temperatures beyond which point vortices form Dirac peaks. 



II. THE TWO-DIMENSIONAL POINT VORTEX GAS 



We consider a multi-species system of point vortices in two-dimensional hydrodynamics described by the Kirchhoff- 
Hamilton equations {33l . l34| : 

dxi dH dyi dH 

dt dyi ' dt dxi ' 

where 7, is the circulation of point vortex i. In an unbounded domain, the Hamiltonian can be written 

ff = E7<7i«(|r<-r i |) > (2) 



i<j 



where 



mr-r'|) = -i-l n |r-r'|, (3) 

is the potential of interaction between two point vortices which satisfies the Poisson equation Au = —S. We note that 
the coordinates x and y of the point vortices are canonically conjugate. In a finite domain, this implies that their 
phase space is bounded (it cannot exceed V N where V is the volume (area) and ./V the number of point vortices). As 
first realized by Onsager this property leads to the occurrence of negative temperature states at high energies. In 
an unbounded domain, the effective size of the system is determined by the conservation of angular momentum (see 
below), so that similar results are obtained. An interesting feature of the dynamics ©-© is that the interaction 
between point vortices determines the velocity of a given vortex while the interaction between material particles (like 
charges or stars) determines the acceleration (force by unit of mass) of a given particle. This is linked to the absence 
of kinetic energy (in the usual sense) in the Hamiltonian . This is another argument why negative temperatures 
are possible for point vortices. For material particles, the temperature is a measure of the average kinetic energy and 
it must be positive (7lj . 

For an isolated system, the energy E = H , the angular momentum L = YliLi 7i r i an d the linear impulse P = 
z x YliLi 7i r are conserved (z is a unit vector normal to the plane of the flow). Of course, the total number N a of 

point vortices of each species is also conserved (the total circulation of species a is T a = N a j a ). The exact velocity of 

a point vortex located in Vi is given by 

Vi = * X fr^ = * X 5>!r7 (|ri " rjl) ' (4) 



where 



V(r)=^ 7lW (|r-r 4 |), (5) 



is the stream function. The stream function is related to the vorticity ui(r) = 7i<J(r — r^) by the Poisson equation 
Aip = —lu. Using Eq. ©, the velocity of a point vortex can be written explicitly 
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We introduce the iV-body distribution P/v(ri, ...,rjv,4) of the system giving the probability density of finding the 
first point vortex in ri, the second point vortex in r 2 ... at time t. The normalization condition is J P/v dri...dr^ = 1. 
We introduce the energy, entropy, angular momentum and impulse functionals 

E[P N ] = J P N H ]Jdri, S[P N ] = -k B J P N \nP N TJdr,, (7) 

N 

L[P N ] = / P N L JTdr,, P[P N ] =Y. PnP Ti dr * ( 8 ) 

i i—1 i 

We define the one- and two-body distribution functions by 

P 1 (a) (r 1 ,t)= f P N (n,...,r N ,t)dr 2 ...dr N , (9) 

P 2 (a&) (n,r 2 ,i) = J P N (v 1 ,...,v N ,t)dr 3 ...dr N . (10) 
The average energy, angular momentum and impulse can then be written 

£=^JV (iV 6 -y J P 2 (ab) (r,r',t) 7a7 ^(|r~r'|)drdr', (11) 

a, b 

L = J2 f N a p[ a \v,t) la r 2 dv, P = -zx^| N a p[ a \v,t) la vdv. (12) 

Remark 1: We note that the angular momentum L = X)i 7i r l i s similar to the moment of inertia for material 
particles. On the other hand, the conservation of linear impulse P = z x J^. 7^ is equivalent to the conservation 
of R = ^iYi which is similar to the center of mass for material particles. The conservation of R for point vortices 
implies that the "center of vorticity" is fixed while, for material particles, the center of mass has a rectilinear motion 
at constant velocity. 

III. THE STATISTICAL EQUILIBRIUM STATE 

We consider a system of point vortices at statistical equilibrium in an unbounded domain. In that case, the 
circulations of the point vortices must have the same sign otherwise they would form dipoles (+, — ) and ballistically 
escape to infinity so that no equilibrium state would be possible (34|. 



A. The YBG hierarchy in the canonical ensemble 

We first consider the statistical equilibrium state of a point vortex gas in the canonical ensemble (CE) in which 
the inverse temperature (3 = l/(ksT) and the angular velocity S7 are prescribed (we set the linear velocity U = 0). 
For Hamiltonian systems with long-range interactions, whose e nerg y is non-additive, the statistical ensembles may 
not be equivalent, even at the thermodynamic limit N — > +00 [35| . Furthermore, for these systems, the canonical 
ensemble has no p hysical sense because it is not possible to define the notion of a thermal bath, especially at negative 
temperatures |72|]. However, it is always possible to consider the canonical ensemble at a mathematical level and study 
the corresponding equilibrium states (the more relevant microcanonical ensemble will be considered in Sec. IIII Cp . 
This study is useful because, at the thermodynamic limit, a canonical equilibrium state is always a microcanonical 
equilibrium state with the corresponding value of the energy. The converse is wrong in case of ensemble inequivalence. 

The canonical iV-body distribution is given by 
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Because of the factor e ^2 L ; the distribution (fl3l) is possibly normahzable only if 

sgn(7)/3fl > 0. (14) 

Using the normalization condition J P/y dri...drjv = 1, the partition function is given by Z(/3,Q) = 
f e-P H -P9 L drx...dr N . The free energy is defined by F(f3,£l) — — jjlnZ((3, f2) and the Massieu function by 
J(/3, f2) = — /3P(/3, Vt) = lnZ(/3, fi). From the canonical distribution (fl~3|) . it is easy to obtain the following expressions 
for the averages and the fluctuations of energy and angular momentum: 

d(/3F) dtfF) 

{H) = —> {L)=2 dWY { ) 

{H 2 )-(H) 2 = k B T 2 ^, (L 2 ) - (L) 2 = -2^y, (16) 

where the independent variables are /3 and /3f2. These expressions show in particular that the specific heat C = dE/dT 
is positive in the canonical ensemble (even at negative temperatures). Finally, we note that the canonical distribution 
(flU)) minimizes the free energy functional F[P/v] = E[Pn] — TS[Pn] + §L[P/v] at fixed normalization. At statistical 
equilibrium, using Eq. (fig]), we get F[P$] = F(/3 7 Q). 
Differentiating Eq. (fT3| with respect to r^, we obtain 

^f = -^ ^^^(Ir.-r.D + ^r^ . (17) 

Taking i = 1, integrating Eq. (|17p on dr2-.-dr]y, and using Eqs. © and (|TU)) . we obtain the first equation of the YBG 
hierarchy in the canonical ensemble 



dp[ a) 

dr 



(r) = -/3£(JV 6 -<U) J P 2 (ah) (r ) r')7a76^(|r-r'|)dr'-/30 7o P 1 (a) (r)r. 



(18) 



B. The virial theorem from the canonical YBG hierarchy 

We introduce the number density n a (r) = N a P^ a \r) of point vortices of species a and the total number density 
n(r) — S(t — rj)) = ^ a n a (r). Similarly, we introduce the vorticity ui a (r) = N a j a P^ a \r) = 7 a n. a (r) of species a 
and the total vorticity uj(r) = 7i<5(r — r-j)) = J2 a w a(r). It is also convenient to define the local "pressure" [l?} by 
the relation 

p{v)=n(v)k B T = k B TY, 1 ^- (19) 

This relation is similar to the local isothermal equation of state of a self-gravitating system in the canonical ensemble. 
For point vortices, the pressure is positive at positive temperatures and negative at negative temperatures. Multiplying 
Eq. (|T5|) by N a and summing on the species, we obtain 

Vp(r) = -J2Na(N b -S ab ) lalb J P 2 (ab) (r,r')|^(|r-r'|)dr'-r! W (r)r. (20) 

ab 

This identity is similar to the- condition of hydrostatic equilibrium for a self-gravitating system in the rotating frame 
(this analogy will become clearer in the mean field approximation, see Eq. (|49[1 ). Taking the scalar product of Eq. 
(|20[) with r and integrating over the whole domain, we obtain 

J \7p-rdr = -V - OL, (21) 

where V is the average value of the virial of the point vortex gas defined by Eq. (|A9[) and L = J cor 2 dr is the average 
angular momentum. Integrating the first term by parts with PV — > at infinity, and using the isothermal equation 
of state (HHJ), we get 

iflL = Nk B T - iv. (22) 
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Finally, using the expression (|A11[) of the virial in an unbounded domain, we obtain the exact virial theorem of point 
vortices at statistical equilibrium in the canonical ensemble 



UlL = Nk B {T-T c ), (23) 



with the exact critical temperature 



k B T c = — — YY 7<7i = — — ( r 2 -Yin = — — (r 2 - r 2 ) , (24) 

8ttN h ' 3 8ttN \ ^ h 8ttN V ; ' y ' 

where T = ^ a N a j a is the total circulation and T2 — J2 a N a lt the total "enstrophy" . For a system of N vortices 
with equal circulation 7, it reduces to 

k B T c = -(N-l)f-. (25) 

We note that the critical temperature is negative. Equation (|23p can be derived directly from the partition function 
(see Appendix IB 1[) . It generalizes the relation previously obtained by [f| [2(| [29| in the single species case and in the 
mean field approximation. Our approach clearly shows that this relation can be interpreted as "the virial theorem of 
point vortices" . Some consequences of the virial theorem are discussed in Sec. IIVEI 



C. The YBG hierarchy in the microcanonical ensemble 

We now consider the statistical equilibrium state of a system of point vortices in the microcanonical ensemble 
(MCE) where the energy E, the angular momentum L and the center of vorticity R are fixed. The microcanonical 
ensemble is the proper description of an isolated Hamiltonian system with long-range interactions [35| . 

The microcanonical W-body distribution is given by 



Ptrin, v N ) = ^rjjHE - H(r u v N ))5 ^ - E 6 ( R - E 



(26) 



Using the normalization condition J Pjy rfri...drjv = 1, the density of states is given by g(E,L) = J S(E — H)S(L — 
J2i 7i r2 )^(R- — Yli Ji r i) dri...drjv- It does not depend on R since this constraint may be absorbed by a shift of origin 
in the integral [6|. We take the origin of the coordinates at the center of vorticity so that R = 0. The entropy is 
defined by S(E, L) = k B ha.g(E, L). Then, the temperature and the angular velocity are given by 

1 = os_ t^^ds 

T dE' 2T dL' [ ' 

Finally, we note that the microcanonical distribution (|26[) maximizes the entropy functional S[Pn] at fixed energy, 
angular momentum, impulse and normalization. At statistical equilibrium, using Eq. (f2l)|). we get S[P^] = S(E,L). 
Differentiating Eq. © with respect to ri and using the microcanonical distribution (|2l)]) . we obtain 



b 



g{E,L) 8E 



g(E, L) 8L (^ P M ^ - yjE~VjdR, (^^(r)) 7- 0») 

This is the first equation of the YBG hierarchy in the microcanonical ensemble 0. This equation involves the 
quantities 

18/ „(ab)\ „„(ab) . dPr ah 



gdE 



(gP^)=0P^ + ^L., (29) 
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The last terms in these expressions are difficult to evaluate. If we neglect these terms, Eq. (|28[) reduces to 



dP{ 



(a) 



Or 



(v) = -l3j2(N b -S ab ) J P { 2 ah \vy) lalb ^{\v-v'\)dv' - pn la p[ a \v)v, (32) 



and we recover the same equation as in the canonical ensemble [compare Eqs. (|32|) and (p~8])] . However, further 
considerations show that the last terms in Eqs. (|2T))) - ([3"l"T) can be neglected only in the N — > +00 limit (see Sec. II VI) . 
Nevertheless, we show below that it is possible to obtain an exact explicit expression of the virial theorem in the 
microcanonical ensemble that is valid for any N . 

D. The virial theorem from the microcanonical YBG hierarchy 

Starting from the hrst equation (|28l) of the YBG hierarchy and using a procedure similar to the one developed in 
Sec. IIIIB1 we obtain 

f Vp ■ rdr = -~(ffV) - ~{2gL) -Jr(gR), (33) 



k B T J r gdE™ ' g dL x * ' g <9R V 

where V is the average value of the virial of the point vortex gas defined by Eq. (|A9[) . L = Jur 2 dr is the angular 
momentum, and R = J ljv dr is the center of vorticity. Expanding the derivatives and using Eq. (|27p , the foregoing 
equation can be rewritten 

dV 

Vp ■ r dr = - V - flL - k B T- Ak B T. (34) 

oE 

Integrating the first term by parts with PV — > at infinity, and using the isothermal equation of state (fl~9|) , we get 

UlL=(N-2)k B T- l -V- l -k B T^. (35) 

Finally, using the expression (|All[) of the virial in an unbounded domain, we obtain the exact virial theorem of point 
vortices at statistical equilibrium in the microcanonical ensemble 

^lL = Nk B {T-T c )-2k B T, (36) 

where we have introduced the critical temperature (|2"4"]) . The same expression can be derived directly from the density 
of states (see Appendix IB 2[) . The result is slightly different than in the canonical ensemble [sec Eq. (l23l) ] because 
of the term —2k B T. Furthermore, for a general potential of interaction u, Eqs. (f2"2"j) and ([33]) differ by an additional 
term —(l/2)k B TdV/dE that may not vanish. However, the canonical and microcanonical expressions coincide in the 
N — > +00 limit. 

IV. THE MEAN FIELD APPROXIMATION 

A. The Boltzmann-Poisson equation from the YBG hierarchy 

We consider the thermodynamic limit N -> +00 with 7 ~ 1/N [uhToT ]. This implies 5~7V, /3~iV, L~l 

and fl ~ 1. Alternatively, one could introduce the dimensionless quantities E* = E/N 2r y 2 : S* = S/N, ft* = Nj 2 (3, 
L,, = L/N^ and 0» = Cl/Nj that can be obtained from simple dimensional analysis (we note that the dimensional 
energy scales like iV 2 ). In this thermodynamic limit, it can be shown that the reduced correlation functions scale like 
Pj ~ {l/Ny^ 1 @, [23. In particular, the two-body correlation function scales like P! 2 ~ 1/N (except in the region 
where |r — r'| is small). Therefore, when N — > +00, the mean field approximation becomes exact: 

N 

P N (r 1> ...,r N )=l[P 1 (r i ). (37) 
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In particular, we have 

P 2 (a6) (r,r') = P^ a \r)P^{r'). (38) 

Furthermore, for N — > +00, the last terms in Eqs. (|29|) -(|31 |) can be neglected, and we can make the approximation 
N a — 1 ~ N a . Therefore, the first equation of the YBG hierarchy in microcanonical and canonical ensembles can be 
written 

VP! (a) (r) = -^"'(^(r) - pj a np( a \r)r, (39) 

where ^(r) = 7iit(|r — r*[)) = J ^ b NbjbPi (r')u(\r — r'|) dr' is the average stream function @. It is related to 
the average vorticity uj(r) = (£\ 7,<5(r - r 4 )) = J2b N blbPi (r) by 

i>(r) = J u(\r-r'\)u(r')dr', (40) 

which is the solution of the Poisson equation 

Aip = -lj, (41) 

with the Gauge condition tp + ^ In r — >• for r — >• +00. The mean velocity of a point vortex is u = z x V?/>. 
Equation (|39p can be integrated to yield the Boltzmann distribution 

W .(r) = AT a7a Pi (a) (r) = r a /e _^ //(r)dr , (42) 

where ip e ff(r) = ip(r) + §r 2 is the relative stream function accounting for the conservation of angular momentum. 
The Boltzmann distribution (142[) is steady in a frame rotating with angular velocity f2. The vorticity profile decreases 
at large distances as 

w„(r) oce-^5 r2 r^, (r -> +00). (43) 

Summing Eq. P2")l on the species and substituting the resulting equation in Eq. (jH]) . we obtain the multi-species 
Boltzmann-Poisson equation 

-A^ = ^r _^ //(r)dr . (44) 

a J 

When /3 > 0, the Boltzmann-Poisson equation (|4"4")) is similar to the one appearing in the theory of electrolytes (like- 
sign vortices "repel" each other) Q and when /3 < 0, it is similar to the one appearing in the statistical mechanics of 
stellar systems (like-sign vortices "attract" each other) 0, [l(| ■ These analogies were noted in @, [H| . We emphasize 
that these results are valid both in microcanonical and canonical ensembles. Furthermore, Eqs. (|42)) and (|44[) remain 
valid in a bounded domain T>, in which case the Boltzmann-Poisson equation must be solved with the boundary 
condition ip e f f = on dT> [f| . In a disk, the angular momentum is conserved, as in an unbounded domain, and we 
get exactly the same equations (with different boundary conditions). If the domain does not possess the rotational 
invariance, the angular momentum is not conserved and £1 — 0. In a bounded domain, we can consider equilibrium 
states of point vortices with positive and negative circulations, since the boundary constrains the phase space to finite 
volumes. 

In a neutral system of point vortices at positive tem per atures, each vortex has the tendency to be surrounded by 
vortices of opposite sign which screen the interaction |36l . [38j . In that case, the system is spatially homogeneous. 
In fact, if we linearize Eq. (|4"4"f for j3 — > (high temperatures), we obtain — (J2 a n alt)P' i l ) = which shows 
that the interaction is exponentially screened on a typical distance A^> = [Q^ a n a r ^)fi\~ 1 l' 1 corresponding to the 
Debye length in plasma physics. At lower temperatures (higher /3), the vortices have the tendency to form pairs or 
dipolcs (+, — ), and the Debye-Hiickel approximation is not valid anymore (see Appendix [E]) . In a neutral system of 
point vortices at negative temperatures, each vortex has the tendency to be surrounded by vortices of the same sign. 
This corresponds to a form of anti-shielding leading to the formation of clusters [H, Hg|. For small /?, we obtain 
^V-' + (J2a n a'~fa)\fi\' l l J = 0, s0 the effect of negative temperatures is to convert the exponentially damped linearized 
solutions discussed above into spatially oscillatory ones [3j. At smaller /3, the system is spatially inhomogeneous and 
consists in large-scale vortices (typically a big dipole). For a neutral system with N/2 vortices of circulation +7 and 
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N/2 vortices of circulation —7, we have uj + = ^4 + e _ ^ 7 ^ e// and w_ = A-e® 1 ^' 11 . If we assume antisymmetry of the 
charge distribution about the center of the box, leading to A + = —A_ = A, we obtain the celebrated sinh-Poisson 
equation Aip — 2Aamh(f3^/ip e f /) [U, [f|. However, there is no fundamental reason why A + = —A- should hold in 
general, so Eq. (|44p is more general. 

Remark 2: Kiessling [39| has shown that there exist an inverse temperature (3q < above which the homogeneous 
solution ip = is the only solution of the Boltzmann-Poisson equation (|44| in the neutral case. In that case, the 
properties of the system are adequately described by the two-body correlation function (see Appendix [D|) . By contrast, 
for j3 < /?o, the mean field equation f|44[) admits spatially inhomogeneous solutions. Such solutions exist only above a 
critical inverse temperature /3 C (see Appendix [E| . 

Remark 3: For point vortices in an unbounded domain, the invariance by rotation of the system leads to a factor 
e~2p~f^r m the vorticity distribution. When the condition (fT4f is satisfied, this factor prevents the dispersion of 
the vortices and confines the system. For rotating self-gravitating systems in an unbounded domain, the equivalent 
factor in the density distribution is e 2^ m ( f2xr ) . This term, which accounts for the effect of the centrifugal force, 
disperses the particles. As a result, there is no statistical equilibrium state for rotating self-gravitating systems in an 
unbounded domain, contrary to the case of point vortices. 

B. The Boltzmann-Poisson equation from the maximum entropy principle 

In the mean field approximation, the entropy, energy, angular momentum and circulation of each species can be 
written 

S= -k B J2j ^ ln ^ dr ' E=^Jutpdr, L = J cur 2 dr, T a =Juj a dr. (45) 

The entropy (|4"51 a) can be obtained by substituting Eq. (|57|) in Eq. (JTJ-b). It can also be obtained from a standard 
combinatorial analysis [HlHI starting from the Boltzmann formula S = kg In W , where W is the number of microstates 
(complexions) corresponding to a given macrostate [73| . Using the Stirling formula for A ;§> 1, we obtain the expression 
(|45l a) of the Boltzmann entropy @, [2l|. The energy (j45]-b) is obtained by substituting Eq. <j38j) in Eq. (HJ) and 
making the approximation A a — 1 ~ N a . 

In the microcanonical ensemble, the statistical equilibrium state is obtained by maximizing the Boltzmann entropy 
5" while conserving the energy E, the angular momentum L, and the total circulation r a of each species 

S(E,L,T a ) = max {S B [u a ] \ E[u] = E, L[lo] = L, T a [u a ] = T a }. (46) 

This maximum entropy principle amounts to determining the most probable state, i.e. the one that is the most 
represented at the microscopic level. Introducing Lagrange multipliers to take the constraints into account, and 
writing the variational problem in the form SS/ks — (3SE — (3^SL — ^2 a a a ST a = 0, we find 0, [2l| that the critical 
points of constrained entropy are given by the mean field Boltzmann distribution (|42[) . 

In the canonical ensemble, the statistical equilibrium state is obtained by maximizing the Boltzmann free energy 
J = S/kB — (3E ~~ P^L (Massieu function) while conserving the total circulation T a of each species 

J(/3, Q, T a ) = max {J B [u a ] \T a [u a \=T a }. (47) 

Writing the variational problem in the form 8 J — a a ST a = 0, we find that the critical points of constrained 
free energy are also given by the mean field Boltzmann distribution (|42[) . Therefore, the series of equilibria (critical 
points) are the same in MCE and CE but the nature of the solutions (maxima, minima, saddle points) may differ 
in each ensemble. When this happens, we speak of ensemble inequivalence (40| . It can be shown that a solution of 
the canonical problem (|47)) is always a solution of the more constrained dual microcanonical problem (|46|) but the 
reciprocal is wrong in case of ensemble inequivalence. 

C. The mean field virial theorem 



Taking the logarithmic derivative of Eq. (|4"2]l . we obtain 

Vw a (r) = -/37 a tJ a (r)VV>e//(r). 



(48) 
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Dividing Eq. (|48|) by 7 Q , summing on the species and introducing the local pressure (|19|) . we get 

\7p + ujVip eff = 0. (49) 

This equation is valid in microcanonical and canonical ensembles. It is similar to the condition of hydrostatic equilib- 
rium for self-gravitating systems in a rotating frame |4lj . Taking the scalar product of Eq. (|49[) with r. integrating 
over the entire domain and integrating by parts the pressure term, we obtain the mean field virial theorem 

2 J pdr-V -QL = 2PV, (50) 

where V is the mean field virial (|A14|) of the point vortex gas, V the area of the domain, and P = 7^7 § V v ' dS the 
average pressure on the boundary of the domain (if the pressure is uniform on the boundary, with value p&, then 
P = pb). Using the isothermal equation of state (|19[) . we can rewrite the virial theorem as 

2Nk B T - V - S1L = 2PV. (51) 

In an unbounded domain, and for axisymmctric flows in a disk, the mean field virial V is given by Eq. (|A17|) . In that 
case, the mean field virial theorem can be rewritten 

PV = Nk B {T-T c )-^QL, (52) 

with the critical temperature 

ksT c = -^. (53) 
For a single species system, the critical temperature reduces to 



k B T c = --L. (54) 

07T 



We sec that the mean field critical temperature ([54]) only differs from the exact critical temperature (|25j) in CE by 
the replacement of N — 1 by N. Of course, they coincide at the thermodynamic limit N —> +00 since the mean field 
approximation becomes exact in that limit. For a neutral system (r = 0), we find that 

k B T c = 0. (55) 

The relation (|52[) is valid in microcanonical and canonical ensembles. It is similar to the virial theorem of a rotating 
2D self-gravitating system. We now consider particular cases of the virial theorem. 

(i) For an axisymmetric flow in a disk with SI = 0, we obtain the equation of state 

PV = Nk B {T-T c ), (56) 

where P = p(R) = n(R)k B T and V = irR 2 . For a single species system, this relation can be directly derived from the 
analytical solution of the 2D Boltzmann-Poisson system (see Appendix IB 3|) . On the other hand, when T = 0, Eq. (f5"6")) 
reduces to PV = Nk B T like for a perfect gas. This relation is of course trivial if the flow is spatially homogeneous. 

(ii) In an unbounded domain, PV — > provided that the pressure decreases sufficiently rapidly with the distance. 
In that case, Eq. ([5^)1 reduces to 

-QL = Nk B (T-T c ). (57) 

This returns the relation obtained in 0, l26l [29j for the single species point vortex gas. However, our derivation 
is valid for a multi-species system of point vortices and does not assume (or use the fact) that the distribution is 
axisymmetric. The mean field virial theorem (|57|) has the same form as the exact virial theorem (|23|) in CE and (|B17j) 
in MCE, but the expression of the critical temperature is slightly different. However, all the expressions coincide at 
the thermodynamic limit N — > +00. 

Remark 4'- We note that the condition L = in MCE does not imply SI = 0. There can be a rotation of the system 
even if the angular momentum is zero. Therefore, we must take the conservation of angular momentum into account 
even when L = 0. 
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D. Summary 

Regrouping the previous results, we find that, in the mean field approximation, the equilibrium stream function is 
determined by the multi-species Boltzmann-Poisson equation (|44[) with appropriate boundary conditions. Then, the 
distribution of each species is determined by the Boltzmann distribution (|4"2"j) . These results arc valid in MCE and 
CE. In CE, the temperature and the angular velocity are prescribed. In MCE, they must be expressed in terms of the 
energy and angular momentum that are the relevant control parameters (conserved quantities) in that case. This is 
done by substituting the Boltzmann distribution (|42j) in the constraints (|45j) . In an unbounded domain, we can use 
the virial theorem (|57|) to relate the angular velocity to the temperature according to 



pL V Pc 

where (3 C = — 8irN/T 2 . This equation replaces the constraint on the angular momentum (|45Ich Then, the temperature 
is determined by the energy (|4"5l b) and the chemical potentials by the circulation of each species f4"5l d). These 
equations generalize the equations given by Lundgren & Pointin [6j for the single species point vortex gas in an 
unbounded domain. The above procedure just determines critical points of entropy in MCE and critical points of 
free energy in CE. We must then study the second variations of these functionals to select maxima of entropy S and 
free energy J. Minima or saddle points must be discarded. When several maxima exist for the same values of the 
constraints, we must distinguish stable (global) and metastable (local) states. Finally, phase transitions and possible 
situations of ensemble inequivalence must be considered. 

E. Consequences of the virial theorem in an unbounded domain 

We consider some consequences of the virial theorem ([57)1 or (|58p in an unbounded domain. In that case, all the 
vortices must have the same sign. To be specific, we assume that their circulation is positive: sgn( / y) > 0. 

1. Microcanonical ensemble 

If L = 0, the vorticity distribution is a Dirac peak at r = containing all the vortices: w(r) = r<5(r), and its energy 
is infinite E — > +oo. This solution exists only at /3 = /3 C (i.e. T = T c ). We now assume L > which is the generic 
case. Because of the inequality (fT4"f and the virial theorem (|55)l . statistical equilibrium states can possibly exist only 
for /3 > j3 c . If j3 = p c , the angular velocity O = 0. In that case, the vorticity profile of species a decreases as r ~ 4Nja / r 
at large distances [see Eq. (|43|) ]. This implies that the angular momentum is infinite (74| so this solution must be 
rejected. Therefore, statistical equilibrium states can possibly exist only for f3 > j3 c (i.e. T > or T < T c ). In the 
single species case, Lundgren & Pointin [(| have shown that the inverse temperature j3 is a function of the augmented 
energy E = 8irE/T 2 + ln(L/r) only. Using the virial identity ([55)1 . we conclude that ilL is also a function of E only. 
Let us consider particular limits. 

(i) If P ->• p+ (i.e. T -> T~), the angular velocity fl ~ -{T 2 /AttL){1 - P/p c ) -> 0". In the single species case, 
Lundgren & Pointin [|| proposed an approximation of the vorticity profile of the form 

W(P ) = _ _2 - e -i (59 ) 

1 + r pj 

where luq is determined by the constraint V = J u dr leading to the relation ttLujo/T 2 + ln(7rLwo/T 2 ) = — 7_e — ln(l — 
PI '/3 C ) where je — 0.557 is the Eulcr constant. The first factor in Eq. ([59jl is the exact solution of the 2D Boltzmann- 
Poisson equation with O = (see Appendix IB 3p and the second factor is a correction for large r that reproduces the 
exact asymptotic behavior (|43[) of the density profile. In this limit, the relation between the inverse temperature and 
the augmented energy is P/P c = 1 — exp[— 1.5772 — E — exp(l + E)]. Therefore, P tends towards p c very rapidly as 
E —> +oo. For p —> /?+, the system develops a Dirac peak at r — containing almost all the vortices, surrounded by 
a residual density profile with a Gaussian tail that ensures the conservation of angular momentum. 

(ii) If p — s- 0* (i.e. T — s- ±oo), the angular velocity ft ~ 2N/PL -t ±oo and Pft — > 2N/L. In that case, the vorticity 
distribution of each species is Gaussian 

/ s N~f a T a _^a r 2 

u a {r) = — e e . (60) 

1TL 
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FIG. 1: Microcanonical parameter space (E,L). The iso-/3 and iso-fiL lines correspond to E = Cst. leading to the family of 
straight lines ln(L/T) = E- 8ttE/T 2 . 



In the single species case, the augmented energy corresponding to j3 = is E(0) = je — hi 2 ~ —0.116. 

(hi) If (3 — >• +oo (i.e. T — >• + ), the angular velocity 17 — > T 2 /(4nL). Since fiip e ff cannot diverge [see Eq. (|32"])]. we 
must have V'e// ~^ 0. This implies ^ ~ ~^r 2 . In that case, the vorticity has a uniform value 



r 2 



(61) 



in a disk of radius tq = (2L/T) 1 / 2 , while it vanishes outside. The augmented energy corresponding to /3 — > +oo is 
E{fi = +oc) = 1/2 - In 2 ~ -0.193. 

Lundgren & Pointin have numerically solved the Boltzmann-Poisson equation (j44]l for a single species point 
vortex gas in an unbounded domain and they have computed all the useful thermodynamic functions [75j . As a 
complement, we have represented the iso-/3 lines and the iso-fiL lines in the microcanonical parameter space (E, L) 
in Fig. [U 



2. Canonical ensemble 



Considering the inequality (I14[) and the virial theorem (|58[) . we come to the following conclusions: 

(i) If Q, = 0, statistical equilibrium states can possibly exist only at T — T c (i.e. f3 = f3 c ). The angular momentum 
is infinite while the energy can take arbitrary values. For a single species system, the vorticity profile is given in 
Appendix IB 31 

(ii) If fi > 0, according to the virial theorem (|58|) . statistical equilibrium states can possibly exist only for T > T c . 
In fact, according to the inequality (fTlj) . when O > the canonical distribution can possibly exist only for T > 
(i.e. /3 > 0) which is a more stringent condition than the one coming from the virial theorem. For ft — > +oo, 
we must have ipeff implying ~ — ^r 2 . In that case, the vorticity has a uniform value 2J7 in a disk of 
radius tq = (r/27rf2) 1 / 2 while it vanishes outside. The angular momentum is L = r 2 /(47rf2) and the energy is 
8irE/T 2 = E(f3 = +oo) - ln(r/47rf2). For f3 -> + , we have L -> +oo and E -> — oo, except if O ^ +oo such 
that f3Q is finite. In that case, the vorticity distribution of each species is Gaussian: ui a = (T a j a (3Q /2ir)e~ la ^ nr I 2 . 
Furthermore, L = 2N/(PQ) and 8nE/T 2 = E(0) - ln(2/7/3f2) for a single species system. 

(hi) If il < 0, according to the virial theorem statistical equilibrium states can possibly exist only for T <T C 
(i.e. j3 c < (3 < 0) which is a more stringent condition than the inequality (fT4"f yielding T < (j3 < 0). If T = T c 
(i.e. /3 = /3 C ), the virial theorem implies L = 0. In that case, the system forms a Dirac peak at r = containing 
all the vortices: ui(r) = TS(r), and the energy E — > +oo. For /3 — ^ 0~, we have L — > +oo and i? — ^ — oo, 
except if fi — > — oo such that /3fi is finite. In that case, the vorticity distribution of each species is Gaussian: 

w a = (r a 7 a /3f2/27r)e- 7 »' 3n '' 2/2 . Furthermore, L = 2N/{fiU) and 8ttE/T 2 = E(0) - ln(2/ 7 /3f2) for a single species 
system. 

We have represented the iso-E lines and the iso-L lines in the canonical parameter space (j3, fi) in Fig. [5] The 
ensembles arc equivalent although the relation between (E, L) and (/3, f2) is not completely trivial. In particular, there 
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FIG. 2: Canonical parameter space (/3, fi). The iso-L lines correspond to $1 = (2N//3L)(1 — f3/f3 c ), and they tend to F 2 /(iirL) 
for ft — > +oo. The iso-E lines correspond to ft = Cst. The intersection between these two curves determines the energy and the 
angular momentum (E,L). 

exist a critical point (O = 0, = /3 C ) in the canonical ensemble corresponding to an infinity of solutions with arbitrary 
energy (— oo < E < +oo) and infinite angular momentum (L = +oo) in the microcanonical ensemble (see Appendix 
IB 3[) . This corresponds to a situation of marginal ensemble equivalence [4(| • 

F. Consequences of the virial theorem for axisymmetric flows in a disk 

We consider some consequences of the virial theorem (|52p for an axisymmetric flow in a disk. 

(i) If the vorticity field is a Dirac peak at r = containing all the point vortices, then L = P = 0, implying T = T c . 

(ii) If = 0, the virial theorem (|52|) reduces to Eq. (|56|) . Since the pressure has the same sign as the temperature, 
an equilibrium state can possibly exist only for j3 > f$ c (i-e. T > or T < T c ). When r = 0, implying T c = 0, the 
virial theorem does not put any constraint on (3. If P = (e.g. if all the vortices form a Dirac peak at r = 0), then 
ft = (3 C (i.e. T = T r ). For the single species case, the vorticity profile is given in Appendix IB 31 

Smith & O'Neil [42( have numerically solved the Boltzmann-Poisson equation (|44|) for a single species point vortex 
gas in a disk by taking the conservation of angular momentum into account. They discussed bifurcations between 
axisymmetric solutions and off-axis solutions, and reported situations of ensemble inequivalence. Caglioti et al. flij 
proved under which conditions the equivalence of canonical and microcanonical ensembles holds. 

Remark 5: In terms of the relative stream function ■0 e //(r) = ip( r ) + § r2 > the Poisson equation can be written 
Atpeff = -lo + 2il. When Q = O c = T/(2ttR 2 ), we have the analytical solution w = r/(?ri? 2 ) and ip = -Tr 2 /(4ti\R 2 ), 
corresponding to ipeff = 0. We can explicitly check that the virial theorem (|52p is satisfied for this solution. In MCE, 
this solution exists only for L = TR 2 /2 and E = — r 2 /(167r). In CE, it exists for ft = £l c at any temperature /3. 

V. CONCLUSION 

In this paper, we have shown that the "Williamson relation" [26| between the angular velocity, the angular momen- 
tum and the temperature of the point vortex gas at statistical equilibrium in an unbounded domain can be interpreted 
as "the virial theorem of point vortices" . Our approach extends the virial identity obtained by Corngold [ItJ which 
is only valid at positive temperatures and applies to a relatively different system. It also extends the calculations 
of Kiessling [29[ that are valid for the single-species point vortex gas in the mean field approximation. We have 
obtained the "exact" expression of the virial theorem and of the critical temperature, valid beyond the mean field 
approximation, for a multi-species system of point vortices. We have discussed the differences between the canonical 
and microcanonical ensembles. We have also derived the mean field virial theorem of an axisymmetric distribution of 
point vortices in a disk although, in that case, the equation is not closed (as it involves the pressure on the boundary). 
Finally, we have discussed some consequences of the virial theorem for the point vortex gas and derived some necessary 
conditions of existence of statistical equilibrium. 
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Appendix A: The virial of point vortices 
1. The exact virial 

Wc define the exact virial of point vortices by 

N 

v = J2w W-Oi). (Al) 

1=1 

This quantity is similar to the virial of a self-gravitating system where the circulation plays the role of the mass and 
the stream function the role of the gravitational potential |4lj . In an unbounded domain, using Eq. ([5]), we obtain 

N „ N _ 

v = EE wi • ^(i* - = EE wi • ^HV^ " rjl) ' (A2) 

where we have used 

S (|r ~ r ' l) = F^ u ' (|r - r ' l} - (A3) 

Interchanging the dummy variables i and j, and adding the resulting expression to the previous one, the foregoing 
expression can be rewritten 



V = \ EE^'h - r>'(| ri - r 3 -|). (A4) 

i=l j^i 

Using the identity u'(|x|) = — l/(27r|x|) coming from Eq. (|3J), we obtain the simple exact result 

N 



iEE™ ( A5 ) 



v 

i=i j^i 

which is independent on the configuration of the system. It only depends on the circulations of the point vortices. 
The sum can be rewritten 

N N 

E E 7i7i = 5>(r - 7<) = r 2 - T 2 , (A6) 

2=1 2=1 

where T = Y^iLi li IS the total circulation and r 2 = Y^iLi if ' IS the total "enstrophy" . For a system of N identical 
point vortices with circulation 7, we get 

V = -N(N-l)-j-. (A7) 

47T 

For a neutral system (r = 0), we have V = r2/(47r). In particular, for a neutral system with AT/ 2 point vortices +7 
and N/2 point vortices —7, we obtain 

V--JL. (AS) 
The ensemble average value of the virial in an unbounded domain is given by 

V = J2"f^N a (N b - S ab ) J P 2 (afc) (r,r') r • ^(|r - r'|) drrfr'. (A9) 

ab 

Using the same trick as before, the foregoing expression can be rewritten 

V = JE^AW-M / P 2 (ab) (r,r')|r-r>'(|r-r'|)drdr'. (A10) 
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Using the identity u'(|x|) = — 1/(2tt|x|) coming from Eq. ((3]), we obtain 



v = J2^N a (N b - s ab ) = -i- (r 2 - r 2 ) , (Ail) 

ab 



where T = J2 a Nala and T 2 = ^2 a N a "f 2 . Of course, this expression can be directly obtained from Eq. (|A5|) . 
The time average of the virial is 

i r N 

(V) = lim -/ dt Y^Mt)- V^(r<(t)). (A12) 
Using the equations of motion it can be rewritten 

i r N 

(V)= lim -/ eft V 7i r i (t).(zxr 1 (t)). (A13) 

Of course, the value of (V) is given by Eq. (|A5|) . However, the point that we want to make here is the following. 
If we integrate Eq. (|A13|) by parts, we get the trivial identity (V) = (V). Therefore, this procedure does not give 
anything, contrary to the similar procedure used for material particles which leads to the usual virial theorem 3l|. In 
the present case, the virial theorem must be derived from the statistical equilibrium state. Indeed, for point vortices, 
the temperature is only defined at equilibrium while for material particles it can be defined at any time as the average 
kinetic energy [3l[ . This implies that the virial theorem of material particles is valid for any steady state, while the 
virial theorem of point vortices is only valid at statistical equilibrium. 

2. The mean field virial 

In the mean field approximation, the virial is 

V = J w(r)r • V^(r) dr. (A14) 

In an unbounded domain, using Eq. (|40p. we obtain 

V = [ w(rWr')r • ^(|r - r'l) drdr' = [ w(rWr> • *~ V V Ar - r'l) drdr' . (A15) 
J Or J |r — r'| 

Interchanging the dummy variables r and r' and adding the resulting expression to the previous one, the foregoing 
expression can be rewritten 

V = i f uj(r)u>(r')\r-r'\u'(\r-r'\)drdr'. (A16) 
Using the identity u'(|x|) = — l/(27r|x|) coming from Eq. ([3]), we obtain 



In the single species case, it reduces to 



In the neutral case (T = 0), we find 



V = -^f. (A18) 



V = 0. (A19) 



In a bounded domain, the virial has not a simple expression in general due to the presence of vortex "images" in 
the potential of interaction. There is, however, an exception. If we consider an axisymmetric vorticity profile in a 
disk, the Gauss theorem, obtained by integrating the Poisson equation (|4T|) , reads 
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where r(r) = JT oj(r')27rr' dr' is the circulation within the disk of radius r. Then, the mean field virial (|A14[) is given 
by 

1 r R dv 1 r R dv 2 r 2 

v =-iL r(r) f *—hh < A21 » 

like in an unbounded domain. 



Appendix B: Alternative derivations of the virial theorem 

1. From the partition function 

In the canonical ensemble, the exact virial theorem (|23p can be directly derived from the partition function Z(/3, f2) 
by making a simple transformation r — > r/y/Q. This is similar to the procedure used by Salzberg & Prager [27], HH 
to derive the exact equation of state of a two-dimensional plasma. From the canonical distribution (|13[) . we find that 
the average value of the angular momentum can be written 

L -—p^r- (B1) 

where 

Z{fi, 0) = J e- iai ln |r '" rj 'e"^ E, ^ dri ... drjVj ( B2) 

is the partition function. Making the change of variables = v^r, (we assume here f2 > 0; the case fl < can be 
treated similarly and leads to the same final result), the partition function can be rewritten 

Ztf, n) = JL e - ^ 7173 1 S'<i 7l7j ln l"l e -l E» dxi ... dxAr . ( B3 ) 

We note that the integral is now independent on ft. Therefore 

i<j 



Substituting this relation in Eq. (|B1[) we obtain 



\tlL = Nk B {T-T c ), (B5) 



with the critical temperature 

ksT c = - 



1 N 

HifEETO (B6) 



These expressions coincide with Eqs. (|2lt)) and (fiM)) . 

For an arbitrary potential of interaction u, the partition function is given by 

Z(fi,(l) = J e-f >H -' 3 * L = fe- p ^^ nMl ' , -' in e-^^dt 1 ...dc N . (B7) 

Making the change of variables Xj = vfh-j, we get 

Z ( J 8 j f2) = JL f e-^^^^-^e-^i^ d^...d^ N . (B8) 



Taking the derivative of this expression with respect to fi, using Eq. (|B1[) . and restoring the original coordinates, we 
obtain 

\toL = Nk B T-\ J dr 1 ...dr^-^^e-^-^ i 5^7<7i|r<-ril« , (|ri-r i |). (B9) 
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In the last term, we recognize the canonical average of the virial (|A4[) . Therefore, Eq. (|B9[) can be rewritten 

^flL = Nk B T - iv. (BIO) 

This coincides with the virial theorem (f2"2")l obtained from the canonical YBG equation. For the potential of interaction 
©, the virial V is given by Eq. (|A5[) . leading to Eq. (|B5I) . 



2. From the density of states 

In the microcanonical ensemble, the virial theorem can be directly derived from the density of states g{E, L) by 
making a simple transformation r — s- sfLx @. Taking the origin at the center of vorticity, so that R = 0, the density 
of states can be written 

g(E,L) = J sIe+^-J^^ 1 ^-^ S^L-J^^S^^ dn...dr N . (Bll) 
Making the change of variables = r j / \/L, we obtain 

g(E,L) = L N ~ 2 J $ f^' + ^^TiTjlnlxi-Xjlj ^1-^7^ S (j2^ x ^ d Xl ...dx N , (B12) 

with 

E' = E+±-lnLj2 li7j , (B13) 

i<j 

where we have used the identity S(Xx) = $( x )- Therefore, the density of states is of the form g(E, L) — L N ~ 2 g(E' , 1) 
and the entropy S(E,L) = k B \ag{E 1 L) satisfies the relation 

S(E,L) = (N -2)k B lnL + S(E',l). (B14) 

Accordingly, 

f)C 1 AC 1 

+ £ 7i7j , (B15) 

i<j 

and 

Using the relations (|27[) , we finally obtain the microcanonical virial theorem 

1 1 N 

-m = {N-2)k B T + — Y J Y,W ( B1? ) 

«=1 jyti 

It slightly differs from the canonical virial theorem (|B5[) - (|B6|) due to the factor (N — 2) instead of A" in front of 
k B T. Our expression (|B17[) also slightly differs from the one obtained by Lundgrcn & Pointin Q which involves 
N — 1 instead of A^ — 2 (apparently, they forgot the contribution arising from the linear impulse). Of course, in the 
thermodynamic limit A" — > +oo, all the expressions agree. 

For an arbitrary potential of interaction u, the density of states is given by 

g{E,L) = J "s(e -X>7^(|r; - r^j S (l - X> r *) * (^ 7m ) dri - drN - ( B18 ) 
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Making the change of variables Xj = r,-/ %/£, we get 



1 {E 1 L)=L N - 2 J 5 ^-^7 ? ; 7j u(VI|x l -x J |)j ,5 ^l-^Ts^j « dx ^- dx N- 



(B19) 



Taking the derivative of this expression with respect to L, using Eq. (|27[) , and restoring the original coordinates, we 
obtain 



m AT n i l 9 

= AT - 2 - - 



2fc B T 



2 8E 



g(E,L) J dr 1 ...dr N -^j-^8 ^E - $^7i7i«(ki ~ r il)j 

V i=l / \i=l / i<i 



(B20) 



In the last term, we recognize the microcanonical average of the virial (|A4|) , Therefore, Eq. (|B20[) can be rewritten 

(B21) 



^-2-I-^-A ( ^, L) v). 



2fc B T 2g(E,L)dE 
Finally, expanding the derivative, and using Eq. (|27[) . we obtain 

iLO = (7V-2)fc B r-iv-ifc B r|^. (B22) 

This coincides with the virial theorem (|35|) obtained from the microcanonical YBG equation. For the potential of 
interaction <j3j) , the virial V is given by Eq. (|A5|) . leading to Eq. (|B17p . 



3. From the 2D Boltzmann-Poisson equation 



For a single-species point vort ex g as with il = 0, the Boltzmann-Poisson equation (|44|) can be solved analytically in 
the axisymmetric case [til. UTI. l26l |29| f37l . In addition to 2D hydrodynamics, thi s analy tical profile appeared in various 



omt vortex ga 

domains such as plasma physics 43|, 0] , astrophysics [9|, I45l - l49j and biology [50l453j . Actually, the 2D Boltzmann- 
Poisson equation is also known as Liouville's equation in differential geometry and the analytical axisymmetric profile 
presented below appears as a special case in a general formula of Liouville [54[ . 

(i) In an unbounded domain, the solutions exist only at the temperature fcsT c = — r7/87r and they are given by 
(see, e.g., (H): 



cj(r) 



w 



r , ( r 2 

-—In hr 

47r V 7rwo 



(B23) 



where ojq is the central density and we have used the Gauge condition -0(r) + ^ lnr — > for i — > +oo. When ujq — > 0, 
the distribution is almost uniform and when ojq — > +oo, it tends to a Dirac peak: w(r) — > r<5(r). Since cj ~ r~ 4 for 
r — > +oo, the angular momentum is infinite (except when w(r) = r<5(r) in which case L = 0). The energy is given by 



E 



r_2 

~8tt 



1 + ln 



(B24) 



and it goes from — oo (as ujq — > 0) to +oo (as uj — > +oo). The caloric curve is a straight line T(E) = T c . The entropy 
is given as a function of the energy by the linear relation 



S = k B \N 



8ttE Arl /7T 
+ AT In — 

r 7 \n 



The free energy J = S — [3E is given by 



J = N 



1 + In 



A r 



(B25) 



(B26) 



18 



and it has the same value for all the solutions, whatever the value of ujq. In MCE, there is only one solution for 
each value of the energy E, and the set of all these solutions has the same temperature T = T c . These solutions are 
stable (global entropy maxima at fixed energy and circulation). In CE, the system exhibits an infinity of solutions 
at T = T c that are parameterized by the central vorticity. There is no solution for other values of the temperature. 
These solutions are metastable. Therefore, if we relax the constraint on the angular momentum (this amounts to 
authorizing L = oo), implying f2 = 0, there is non-uniqueness in the canonical ensemble while there is uniqueness 
in the microcanonical ensemble. This is a situation of marginal ensemble equivalence [40| . The entropy S(E) is a 
straight line, so it is both concave and convex. The inverse specific heat d 2 S/dE 2 = d/3/dE vanishes. The line of 
microcanonical solutions E £] — oo, +00 [ is mapped to a single point (3 = f3 c in the canonical ensemble. For a recent 
example of marginal ensemble inequivalence in the context of fluid mechanics, see (55j . 

(ii) In a bounded domain, the solutions exist only for T > or T < T c = —Y^/%-k (i.e. (3 > (3 C = Sir/V"/) and 
they are given by (see, e.g., [Hj]): 

^) = 4o , i^ -r ^ (B27) 



rR 2 1 - T c /T 



1 



T-T c 



ip(r) 



2k B T 

7 



In 



T c T c rr 

1 + — I - 

T T \R 



(B28) 



where we have used the Gauge condition ip(R) = 0. From Eq. (|B27[) . writing P = p(R) = uj(R)kBT '/ '7 and V = irR 2 , 
we obtain 



PV = Nk B {T-T c ), (B29) 

and we recover the equation of state The central vorticity is given by ujq = (r/7ri? 2 )(l — T c /T) _1 and the 

angular momentum by L = TR 2 (T/T C - 1)[-(T/T C ) ln(l - T c /T) - 1]. When T > 0, the vorticity increases with the 
distance and the vortices tend to accumulate at the boundary (the interaction is "repulsive" like between like-sign 
charges in electrolytes). When T < 0, the vorticity decreases with the distance and the vortices tend to cluster at the 
center of the domain (the interaction is "attractive" like between stars in galaxies). When T — > ±00 (i.e. (3 = 0), the 
distribution is uniform. When T = (i.e. (3 — > +00), the distribution is a Dirac peak at r = R. When T = T c (i.e . 
(3 = f3 c ), the distribution is a Dirac peak at r = 0. Some vorticity profiles are represented in Figs. 4 and 5 of [371 ]. 
The energy is given by 



8tt 



T-T c 



(B30) 



and it goes from (for T = 0) to +00 (as T — > T c ). The transition between positive and negative temperatures 
(corresponding to T —> ±00, i.e. (3 = 0) takes place at the energy E = r 2 /(167r). The entropy is given by 



S = ks 



I y) +2N + n(2^-1) ln(l-^ 



(B31) 



Eliminating T between Eqs. JB30|) and (|B31[) . we get S(E). The caloric curve T(E) and the entropy curve S(E) are 
plotted in Figs. 6 and 7 of [331 - Finally, the free energy J = S — f3E is given as a function of the temperature by 



J 



(yW + ^(g-l)ln(l-| 



(B32) 



In MCE, the solutions exist for any E and they have an inverse temperature (3 > f3 c . In CE, the solutions exist for 
a- n y f3 _ (3 C . When the angular momentum constraint is relaxed (£7 = 0), the ensembles are equivalent. 



Appendix C: Virial theorem in a bounded domain for purely logarithmic interactions 

We consider a system of point vortices in a bounded domain of area V. We give here a non-rigorous derivation 
of the virial theorem by ignoring the contribution of "vortex images" due to boundaries. In other words, we use 
the potential of interaction ([3]) as in an unbounded domain. This approximation has also been made by Edwards 
& Taylor [38| . It is difficult to say to which extent this approximation is accurate. However, one motivation of the 
following calculations is to help interpreting the curious notion of "pressure of point vortices" that we have introduced 
heuristically in Sec. 1111 Bl On the other hand, the following results are exact for particles interacting by a purely 
logarithmic potential u(\r — r'|) = — ^- In |r — r'| in 2D, like plasmas and self-gravitating systems. 
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1. From the YBG hierarchy 



If we ignore the contribution of vortex images, the first equation of the canonical YBG hierarchy is given by Eq. 
(fT5)l . If we define the local pressure by Eq. (jT9")l . we obtain Eq. (HH). Integrating the first term by parts, introducing 
the average pressure 



on the boundary of the domain, we get 



2V 



PV = Nk B T 



pr ■ dS, 



1 1 

2 V -2 nL - 



(CI) 



(C2) 



This is the equivalent of the usual virial theorem for material particles [311 ]. However, for point vortices, it cannot 
be directly derived from the equations of motion (see the discussion at the end of Section |A lj) . Using the expression 
15[) of the virial (valid when vortex images are ignored) we obtain 



PV = Nk B (T-T c ) 



1 



-ClL, 



where T c is the critical temperature 



1 N 



TV 



i=l j=jti 



For a neutral system 



8nN 



k B T c = 



8ttN 



r 2 -r 2 



r 2 



8?riV 

For a system consisting of N + vortices of circulation +7 and 7V_ vortices of circulation —7, we get 

- 2 r (iV+-iV_) 2 ' 



kBTc = h 



1 - 



N 



(C3) 



(C4) 



(C5) 



(C6) 



In particular, for a neutral system with N + = N/2 vortices +7 and iV_ 
is 

k B T c 



7_ 
87r' 



N/2 vortices —7 , the critical temperature 

(C7) 



On the other hand, for a single species system with N vortices 7, we get 

k B T c = -(N -1)1-. 

If the domain has not the rotational symmetry, then O = 0. In that case, the virial theorem reduces to 

PV = Nk B (T-T c ). 



(C8) 



(C9) 



Therefore, the virial theorem determines the exact equation of state of a 2D system with logarithmic interactions. 
Since PV/Nk B T > 0, we conclude that a statistical equilibrium state can possibly exist only if (3/f3 c = T c /T < 1. In 
particular, for a single species system for which k B T c = — (N — 1)^ 2 /8tt is negative, we must have (3 > (3 C (T > or 
T < T c ). On the other hand, for a neutral system for which k B T c = T2/8ttN is positive, we must have (3 < (3 C (T < 
or T > T c ). In the mean field approximation, we recover the results of Sec. IIV CI 

The preceding results are valid in the canonical ensemble. In the microcanonical ensemble, starting from Eq. 
and using the same procedure as before, we obtain the virial identity 



PV = (N -a)k B T- 



-V 



2 kBT dE 



-CIL, 



(CIO) 



with a = 1 in a disk [76[ and a = in a domain that does not have the rotational symmetry (CIL — 0). Using the 
expression (|A5[) of the virial (valid when vortex images are ignored), we obtain 



PV = Nk B (T - T c ) - -CIL - ak B T, 



(Cll) 



where T c is the critical temperature (|C4 
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2. From the partition function 

In the canonical ensemble, the pressure is defined by 

d_F_l_dJnZ_ 

F dv ~ p ov ' (C12) 

where 

Z(fi,V)= I ' e&i:^™^**-^ dr,...dv N , (C13) 



is the partition function (ignoring the contribution of vortex images). We have assumed that the system has not the 
rotational symmetry so that = 0. If we make the change of variables Xj = Ti/y/V, the partition function becomes 

Z((3,V) = V N e^ lnV ^^ f e^<i™> hllXi - x > l dx 1 ...dx N , (C14) 

where the integral is independent on V. From this expression, we get 

^T=y + «^ (C15) 

i<j 

Substituting this relation in Eq. (|C12[) . we obtain the equation of state 

PV = Nk B {T-T c ) 1 (C16) 

with the critical temperature (|C4[) . This result extends to the multi-species case the calculations of Salzberg & Prager 
0|2|. 

For an arbitrary potential of interaction u, the partition function is given by 

Z(/3,y) = | e -^dr 1 ...dr w = y' e ^^<, 7 ' 7jU(|r '^ l) dr 1 ...dr w . (C17) 

Making the change of variables Xj = r»/ VV, we get 

Z(J3,V) = V N [ e-^<^ i7iu(VVlxi -* il) d Xl ...dx N . (C18) 

Taking the derivative of this expression with respect to V, using Eq. (IC12I) . and restoring the original coordinates, 
we obtain 

PV = Nk B T-^J dv^.dvN-^—e-^Y.lil^i-^M'^i-^)- (C19) 

In the last term, we recognize the canonical average of the virial (|A4[) . Therefore, Eq. (|C19[) can be rewritten 

PV = Nk B T -iv. (C20) 

Using the expression (|A5|) of the virial (valid when vortex images are ignored) we recover Eq. (|C16I) . We note that 
Eq. (|C20j) coincides with the virial theorem (|C2j) obtained from the YBG equation. Therefore, we conclude that the 
pressure P defined by Eq. (|C1|) with the isothermal equation of state (fli?]) can be identified with the thermodynamical 
pressure. This was not obvious a priori. Therefore, even if the approach developed in this section is not rigorous 
(since it ignores vortex images) , it shows a "consistency" between the different definitions of pressure that we have 
given in the paper. 
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3. From the density of states 

In the microcanonical ensemble, the pressure and the temperature are defined by 

dS i_di 

where S(E, V) = k B \ng(E, V) is the entropy and g(E, V) the density of states 

g(E,V)= [ s(E+^Y,Wi ^\r i -r j \)dr 1 ...dr N . (C22) 

Wc have assumed that the system has not the rotational symmetry so there is no constraint on the angular momentum. 
With the change of variables x, = r j/ VV, the density of states can be rewritten 



g(E,V) = V N f s(e' + ^-^7 i 7 J ln|x j -x i |)dx 1 ...dx^, (C23) 
J \ *" i<j ' 

where 

E' = E+^-\nVj2%lr (C24) 



47T 



Therefore, the density of states is of the form 

g{E,V)=V N g(E',l), (C25) 



and the entropy is of the form 
From this expression, we obtain 

and 

Using Eq. ([C2T) we get 



S{E,V) = Nk B lnV + S(E',l). (C26) 



l^^ + lf^-Ly^, (C27) 



i -§<*■•>■ < C28) 



PV = Nk B (T-T c ), (C29) 
whh the critical temperature (|C4[) . This result extends to the multi-species case the calculations of Edwards & Taylor 

For an arbitrary potential of interaction it, the density of states is given by 



g(E,V)= f 5{E-H)dv 1 ...dr N = f s(e - ^ 7 , 7jU (|r 4 - ry|)J dn...dv N . 



(C30) 

Making the change of variables Xj = ri/VV , we get 

g(E,V) = V N J 6\^E-^iaXVV\r i -r j \)\dx 1 ...d XN . (C31) 

Taking the derivative of this expression with respect to V, using Eq. (|C21[) . and restoring the original coordinates, 
we obtain 



PV AT 1 1 d 

N 



k B T 2g{E,V)dE 



g(E,V)J dr 1 ..Ar N -^^S(E-H)Y / ^ lj \r l -r J \u / (\r l -r J \). (C32) 
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In the last term, we recognize the microcanonical average of the virial (|A4|) . Therefore, Eq. (|C32[) can be rewritten 

S = N -\^v)^ E ^- (C33) 

Expanding the derivative, and using Eq. (|C21j) . we obtain 

PV = Nk B T - l -V -\k B T^. (C34) 

Using the expression (|A5|) of the virial (valid when vortex images are ignored) we recover Eq. (|C29|) . We note that Eq. 
(|C34[) coincides with the virial theorem (|C10[) obtained from the YBG equation. This allows us to identify expression 
(|C1|) with the thermodynamical pressure. 

Finally, we note that the equations of state (|C16|) and (|C29|) have the same form in the microcanonical and canonical 
ensembles (recall, however, that in the microcanonical ensemble T = T(E) is a function of the energy while in the 
canonical ensemble T is prescribed). This equivalence is not expected to be always true because, for an arbitrary 
potential of interaction u, the expressions (|C20[) and (|C34|) differ by a factor — (1 /2)k B TdV / dE that may not vanish. 
However, in the limit ./V — > +00, this term is negligible and the equations of state always coincide in the microcanonical 
and canonical ensembles. 



Appendix D: The case of a neutral spatially homogeneous system 

We consider here a neutral system made of N/2 point vortices of circulation +7 and N/2 point vortices of circulation 
—7. We assume that the system is spatially homogeneous. The density of each species is therefore n + = n_ = n/2 
where n = N/V is the total density. As shown in Appendix [C] the virial theorem leads to the equation of state 

PV = Nk B (T-T c ), (Dl) 

with 

k B T c = (D2) 

<37T 

Let us first consider positive temperatures. In that case, the system is always spatially homogeneous. The exact 
equation of state (|D1|) - (|D2|) was originally derived by Salzbcrg & Prager [27|, [28| and May (5(| for a spatially homoge- 
neous 2D neutral plasma. According to Eq. (|D1[) . the pressure is negative for T < T c , which is not possible at positive 
temperatures. Actually, it can be shown that the partition function is convergent only for T > T* = j 2 /(4tt) = 2T C 
[57H621 ] (see Appendix [EJ so the equation of state (|D1[) is not valid for T < . For T < , there is no equilibrium 
state and the system "collapses". This leads to the formation of N/2 non-interacting singular pairs (+, — ) of opposite 
charge, corresponding atoms (+e, — e) in plasma physics or to dipoles (+7,-7) in 2D hydrodynamics. For T = T*, 
we have PV = (l/2)NkBT lf . This is the equation of state of an ideal gas of N/2 non-interacting singular pairs (+, — ). 

We write the two-body distribution function as (|r — r'|) = Pq + P^^flr — r'|) where Pq = 1/V is the 

spatially homogeneous one-body distribution and P^ ab \\r ~ r'Q is the two-body correlation function. Due to the 

symmetry of charges +7 and -7, we have P2 ++) = 5 = -P^ + ^ = -P2~ +) = P oH\ r ~ r 'D- Usin g tlie 

second equation of the YBG hierarchy and implementing a Dcbye- rliickel approximation which consists in neglecting 
three-body correlations, it can be shown that the two-body correlation function Hdh{\y ~ r '|) is the solution of the 
differential equation |36|, [H, [64| : 

Ah DH - k 2 D h DH = /3 7 2 ^(x), (D3) 

where ku = (/Jn^ 2 ) 1 / 2 is the Debye wavenumber. The Dcbye- Hiickcl approximation is valid at sufficiently high 
temperatures T (small /3). In particular, it does not signal the collapse of pairs (+, — ) at T < T*. The solution of Eq. 
(ID31) is 

n(2ir) 2 h DH (k) = D , h DH (x) = -^-K {k D x). (D4) 

This result can also be obtained from a diagrammatic expansion [dH [d5| in terms of the plasma parameter e = 
/3j 2 /2ir <C 1. A more precise expression of the correlation function (obtained from a straightforward renormalization 



23 



of the short-range behavior of h) is given by h(x) = 1 — e eK oi,k D x) ^^ ^^ According to Eq. (|D4[) . a vortex of a given 
sign is surrounded by vortices of opposite sign which screen its interaction on a distance of the order of the Debye 
length ~ k^ 1 . 

Let us now consider the case of negative temperatures. In that case, the system is spatially homogeneous above an 
inverse temperature (3q < (see the Remark at the end of Sect. IIV A[) . The exact equation of state (|D1|) - (|D2p is still 
valid [38| and we note that the pressure is negative (which is allowed at negative temperatures). The Debye- Hiickel 
approximation can also be extended to negative temperatures in the range Pq < P < 0. The two-body correlation 
function hjjH (|r — r'|) is the solution of the differential equation (36j : 

Ah DH + k 2 h DH = /3 7 2 ^(x), (D5) 

where kj = (—Pnj 2 ) 1 / 2 is the equivalent of the Jeans wavenumber in astrophysics (66j . The solution of this equation 
is 

k 2 Bi 2 
n{2n) 2 h DH {k) = fc2 _ J fc , , h DH (x) = ^Y (kjx). (D6) 

The oscillatory behavior of the correlation function hjjH{x), or the fact that the denominator of }idh(x) is negative 
for k < kj, is related to the clustering of point vortices at negative temperatures (36J. This is consistent with the 
results of Edwards & Taylor [38} obtained in a different manner. For smaller values of the inverse temperature P, the 
system is spatially inhomogeneous (forming a large-scale dipole) and we must use the mean field results of Sect. IIVI 
that are exact at the thermodynamic limit. For P < /3* = — IGtt/Nj 2 , there is no equilibrium state anymore and the 
system forms a singular macroscopic dipole (f +,f-) (see Appendix |E|) . 
We now consider the virial equation [31] : 

1 f + °° 

PV = Nk B T - -n 2 j 2 V / g(x)xu'(x)2irxdx, (D7) 
4 Jo 

where u(x) is the potential of interaction ([3]). This equation can be directly obtained from Eqs. (|C2|I and (|A10|) . 
and it leads to Eq. (|D1|) . If we replace the correlation function g{x) = 1 + h(x) by the Debye-Hiickel expression 
hDH(x) given by Eqs. (|D4|I and (|D6|) . and carry out the integration, we obtain Eq. (jDip . Therefore, the virial 
equation calculated with the Debye-Hiickel correlation function returns the exact equation of state of a 2D plasma or 
a 2D point vortex gas. However, this is essentially a coincidence since Eq. (|D7[) is generally valid for the correlation 
function g(x), not for Iidh{x)- 



Appendix E: Existence of statistical equilibrium states 

We first consider a system of N point vortices of equal circulation 7 in a bounded domain. At positive temperatures 
P > 0, like-sign vortices "repel" each other and accumulate on the boundary of the domain. As the inverse temperature 
increases (or as the energy decreases), the vortices are more an more concentrated on the boundary. This regime of 
low energies corresponds to a repulsive interaction between vortices, like between electric charges of the same sign in 
electrolytes @- At negative temperatures P < 0, like-sign vortices "attract" each other and form clusters. As the 
inverse temperature decreases (or as the energy increases), the cluster is more and more condensed. This regime of 
high energies corresponds to an attractive interaction between vortices, like between stars in a galaxy |<J[l(|. It can 
be shown [Til IT2I |62| that statistical equilibrium states exist (the partition function is finite) if, and only, if 

< E1 > 

When P +00, all the vortices concentrate on the boundary. When ft fti , all the vortices collapse and form 
a Dirac peak w(r) = r<5(r). In that case, E — > +00. In an unbounded domain, the picture is slightly different (see 
Sec. IIV El for more detailed results in the mean field limit). At positive temperatures p > 0, like-sign vortices "repel" 
each other but the spreading of the vortices is prevented by the conservation of angular momentum in MCE or by 
the confining potential tpconf = m CE. When P — > +00, the vorticity profile is a step function (in the mean field 

limit). At negative temperatures Pi < P < 0, like-sign vortices "attract" each other and form clusters. As the 
inverse temperature decreases, the cluster is more and more condensed. When P — > Pi', the vorticity profile is a 
Dirac peak in CE and a Dirac peak surrounded by a residual vorticity profile with a Gaussian tail (in the mean field 
limit) which allows to satisfy the conservation of angular momentum in MCE. In that case, E — > +00. 
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We now consider a neutral system made of N/2 point vortices of circulation +7 and N/2 point vortices of circulation 
—7 in a bounded domain. At positive temperatures B > 0, the interaction between like-sign vortices is "repulsive" 
and the interaction between opposite-sign vortices is "attractive" . This is similar to the case of electric charges in 
electrolytes Q. At high temperatures (low /3), each point vortex is surrounded by a cloud of opposite-sign vortices 
which screens the interaction (3(| HH . This is similar to the Debye shielding in plasma physics @ . In that case, 
the system is spatially homogeneous and fully "ionized". At low temperatures (high f3), the vortices tend to form 
pairs of microscopic dipoles (+, — ) similar to "atoms" (+e, — e) in plasma physics. The system consists in a spatially 
homogeneous gas of non-interacting microscopic dipoles. This corresponds to the regime of low energies. At negative 
temperatures f3 < 0, the interaction between like-sign vortices is "attractive" and the interaction between opposite-sign 
vortices is "repulsive" . At small /3, the system is spatially homogeneous in average but the vortices tend to cluster. 
This corresponds to a form of anti-shielding in which each point vortex is surrounded by a cloud of similar vortices 
[36|, [H[. This is analogous to the Jeans gravitational clustering in astrophysics [6(|. At lower B, the vortices form 
two (or more) macroscopic clusters (dipoles, tripoles,...) made of vortices of the same sign p], Q. In that case, the 
system is spatially inhomogeneous and the structures can be described within the mean field theory which is exact. 
This corresponds to the regime of high energies. It can be shown [62L IoTI] that statistical equilibrium states exist (the 
partition function is finite) if, and only, if 

4-» = -^</ J <4 + '=^. (E2) 

When B — > /3* , the system forms a singular macroscopic dipole (y+> %~ ) with one Dirac peak +^7 containing 
all the positive vortices and one Dirac peak —-^7 containing all the negative vortices. In that case, E — > +00. 

The critical temperature /3* ^ can be understood by considering the statistical mechanics of a single species system 
of N/2 vortices (since each peak contains N/2 vortices). It can be therefore deduced from Eq. (|E1|) by making the 

substitution N — s- N/2. When B —t /3* , the system forms N/2 singular microscopic dipoles (+, — ) where the vortices 
of each pair have fallen on each other. This leads to N/2 Dirac peaks made of two vortices of opposite circulation. In 

that case, E — > —00. The critical temperature /3* can be understood by considering the statistical mechanics of only 
two vortices of opposite sign. It can therefore be deduced from Eq. (|E1[) by taking N = 2 (since only two vortices 
are involved in the dipole) and reversing the sign (since the attraction of like-sign vortices at negative temperatures 
corresponds to the attraction of opposite-sign vortices at positive temperatures). 

For a multi-species system of point vortices in a bounded domain, the conditions of existence of the statistical 
equilibrium state, and the corresponding expressions of the critical temperature, are more difficult to obtain. Some 
interesting mathematical results have been obtained by Ohtsuka et al. [68| in the mean field approximation. Explicit 
expressions of the critical temperature have also been obtained by Sopik et al. when the system consists of two 
types of like-sign vortices (this study is performed in the context of self-gravitating systems, but it can be readily 
applied to point vortices at negative temperatures). 
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